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Abstract. The Arnowitt-Deser-Misner (ADM) equations are deeply in- 
i \ tertwined with discrete spectral resolutions of an elliptic operator of Laplace 

type associated with the spacelike hypersurfaces which foliate the space-time 
■^j- ■ manifold, and the non-linearities of the four-dimensional hyperbolic theory 

\q \ are mapped into the potential term occurring in this operator. The ADM 

O . equations are here re-expressed as a coupled first-order system for the in- 

£N| ' duced metric and the trace-free part of the extrinsic-curvature tensor, and 

7— I ' their formulation in terms of integral equations is studied. 

O 

o ■ 

The canonical formulation of general relativity relies on the assumption 
i , that space-time (M, g) is topologically ExR and can be foliated by a family 

H ■ of spacelike hypersurfaces St, all diffeomorphic to the three-manifold S. The 

Zv' space-time metric g is then locally cast in the form 

> : 

g = -(N 2 - NiN r )dt ®dt + Ni(dx l <g> dt + dt dx l ) + h^dx 1 % dx 3 , (1) 

5h ' 

where N is the lapse function and N l are components of the shift vector 
of the foliation [1-3]. The induced metric hij on Sf and the associated 
extrinsic-curvature tensor Kij turn out to obey the first-order equations 
[1-3] 

Bh- 

-^=-2NK ij + N ilj + N j[i , (2) 



dK, 



dt 



N W + N 



+ 



N m K ij{m + N^K jm + N^K b 



(3) 



where the stroke denotes covariant differentiation with respect to the induced 
connection on Sj. Moreover, the constraint equations hold. For Einstein 



1 



theory in vacuum they read 

K u ]i -(tTK) v ^0, (4) 

(3) i? - KijK ij + (trK) 2 « 0, (5) 

where ~ is the weak-equality symbol introduced by Dirac to denote equa- 
tions which only hold on the constraint surface [3,4]. From now on, we shall 
impose the vacuum Einstein equations in four dimensions, so that in the 
following equations ^ R%j can be set to zero. 

Equation (3) is quasi-linear in that it contains terms quadratic in the 

extrinsic-curvature tensor, i.e. N (^K^^tK) — 2K im KJ n ^j . We are now aim- 
ing to obtain from eqs. (2) and (3) another set of quasi-linear first-order 
equations. For this purpose, we contract Eq. (2) with K lJ and Eq. (3) with 
h % i . This leads to (hereafter —A = — |J* is the Laplacian on S t ) 

Rij d^ = _ 2NK . jK ij + 2 tfVJV ib -, (6) 



- AN + N 



^R+(tiK) 2 -2K im K i 



+ 



N m {trK)\ m + 2N^Kl 



It is now possible to obtain an evolution equation for (trfsT), because 



a »*. Kij+h <i dK >, 



at 



dt 



(7) 



(8) 



The second term on the right-hand side of Eq. (8) is given by Eq. (7), 
whereas the first term is obtained after using the identity 



(h ij h jt ) 



which implies 



and hence 



dt 
dtii _ 

80 

dt Kl3 



d , 
dt 1 



-h ip h jl 



o, 



9h p i 

dt ' 



dt 



(9) 
(10) 

(ii) 
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By virtue of Eqs. (6)-(8) and (11), and imposing the constraint equation 
(5), we get 

^(trK) = - (A - K^) N + N m {trK\ m , (12) 
which can also be cast in the form (here V m = i m ) 

{m ~ Nm ^ m ) {trK) = (" A +Ki i Rij ) N - (13) 

Note that, without imposing the vacuum Einstein equations and the con- 
straint (5), the right-hand side of Eq. (13) would read instead 

(- A + (tiK) 2 - ^Rijh ij ) N. 

Equation (13) is a first non-trivial result because it tells us that, given the 
operator of Laplace type [5] 

V = - A +K ij K i i, (14) 

the right-hand side of Eq. (13) is determined by a discrete spectral resolution 
of V . By this one means a complete orthonormal set of eigenfunctions 
belonging to the eigenvalue A, so that, for each fixed value of t, the lapse 
can be expanded in the form 

N(x,t) = J2Cxf P x (x,t), (15) 

A 

with Fourier coefficients C\ given by the scalar product C\ = (^f^,N^j. This 

point is simple but non-trivial: for each fixed value of t, the restriction of the 
lapse to Yi t becomes a function on S t only, and hence the Fourier coefficients 
in (15) read 

Cx= f (f p x )*NVhd 3 x = Cx, t , (16) 
J s t 

where the star denotes complex conjugation. In the application to the initial- 
value problem, one studies the lapse at different values of t, and hence it 
is better to write its expansion in the form (15), with the time parameter 
explicitly included. The integration measure for C\, however, remains the 
invariant integration measure on S t . 

The existence of discrete spectral resolutions of V is guaranteed if E$ is 
a compact Riemannian manifold without boundary [5], which is what we 
assume hereafter. Thus, Eq. (13) can be re-expressed in the form 

- N m V m ^j trK(x, t) = ^ XC\ t f p x (x, t). (17) 



3 



On denoting by L the operator 




N m V. 



(18) 



and writing G(x, t; y, t') for the Green function of L, Eq. (17) can be there- 
fore solved for (trif ) in the form 



where dfj,(y, t') is the invariant integration measure on T, t i x R. Note that the 
eigenfunctions f% are different at different times and hence, to compute them 
at all times, one has to integrate the equations of motion so that the metric 
and the extrinsic-curvature tensor are known at all times. Thus, Eq. (19) 
does not reduce the amount of entanglement of the original equations (2) 
and (3). A relevant particular case is obtained on considering the canonical 
form of the foliation [6], for which the shift vector vanishes. The operator L 
reduces then to J^. 

A cosmological example shows how the right-hand side of Eq. (17) 
can be worked out explicitly in some cases. We are here concerned with 
closed Friedmann-Lemaitre-Robertson- Walker models, for which the three- 
manifold S reduces to a three-sphere of radius a. It is indeed well known 
that the space of functions on the three-sphere can be decomposed by us- 
ing the invariant subspaces corresponding to irreducible representations of 
0(4), the orthogonal group in four dimensions. These invariant subspaces 
are spanned by the hyperspherical harmonics [7], i.e. generalizations to the 
three-sphere of the familiar spherical harmonics. The scalar hyperspherical 
harmonics Q^ n \x-, f) form a basis spanning the invariant subspace labeled 
by the integer n > 1 corresponding to the n-th scalar representation of 0(4). 
The label (n) of refers both to the order n and to other labels denoting 

the different elements spanning the subspace. The number of elements Q/ n -* 
is determined by the dimension of the corresponding 0(4) representation. 
On using the general formula for 0(2k) representations, the dimension of the 
irreducible scalar representations is found to be d,Q(n) = n 2 . The Laplacian 
on a unit three-sphere when acting on 0/ n ) gives 



Any arbitrary function F on the three-sphere can be expanded in terms of 
these hyperspherical harmonics as they form a complete set: 




AQ>) = (n 2 - 1 



)Q (n) 



n 



(20) 



(21) 



(n) 
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where the sum (n) runs over n = 1 to oo and over the n 2 elements in each 
invariant subspace, and the are constant coefficients. The lapse function 
on the three-sphere is therefore expanded according to (cf. (15)) 

N = ^c in) Q {n \ (22) 

(n) 

and bearing in mind that K^K 13 = 3 on a unit three-sphere we find, on a 
three-sphere of radius a, 

VN = Y,c {n) \ n Q {n \ (23) 

(n) 



where (here a = a(t)) 



X n = { ^±^ n = l,2,... . (24) 



We have therefore found that all nonlinearities of the original equations give 
rise to the potential term in the operator V defined in Eq. (14), which is 
a positive-definite operator of Laplace type (unlike the Laplacian, which is 
bounded from below but has also a zero eigenvalue when n = 1). 

The interest in the evolution equation for the trace of the extrinsic- 
curvature tensor is motivated by a careful analysis of the variational problem 
in general relativity. More precisely, the "cosmological action" is sometimes 
considered, which is the form of the action functional I appropriate for the 

case when (irK) and the conformal three- metric = h~s hij are fixed on 
the boundary. One then finds in c = 1 units [8] 

1 - lib /„ " ,R ^ d ' x + lib U tlK ^ h d? ' x ' < 25) 

and this formula finds important applications also to the quantum cosmology 
of a closed universe, giving rise to the if-representation of the wave function 
of the universe [9]. 

The program we have outlined in our paper shows an intriguing link 
between the Cauchy problem in general relativity on the one hand [10-13], 
and the spectral theory of elliptic operators of Laplace type on the other 
hand [5], which is obtained by exploiting the nonlinear form of the ADM 
equations and the nonpolynomial form of the Hamiltonian constraint (5). 
On denoting by Oij the trace- free part of the extrinsic-curvature tensor: 



<Tij = Kij - -hijitvK) (26) 
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we therefore obtain a set of equations, equivalent to (2) and (3), in the form 
(here L _1 denotes the inverse of the operator L defined in (18)) 

trK = L- 1 VN, (27) 

= --NhijitrK) - 2Na ij + 2N m , (28) 

d ^L = n t] N-^N m (tTK) + N m V m a lJ + 2N^[a j)m + h^^ (29) 
where the operator of Laplace type Qij is given by 

Ojj = -VjVi - -hijV + -aij(tvK) 

+ ^h ij (tiK) 2 + ^R ij -2a im a J m . (30) 

Hereafter, to be more rigorous, we assume that the left and right inverses 
of L exist and are the same, and that left and right inverses coincide for all 
operators we have to consider. 

It now appears desirable to understand to which extent the differential- 
operator viewpoint plays a role in solving the coupled system (27)-(30). 
For this purpose, we consider a relevant particular case, i.e. the canonical 
foliation studied in Ref. [6], for which the shift vector vanishes. The operator 
L reduces then to and, on defining the operators 

Q = §1 + \ N ^ K = L + \NL- l VN, (31) 

ON N 
S=---KK = L--L-^N, (32) 

Aij = .Vj-Vi + (3) %, (33) 

one finds the equations 



2 



2 



5+-(L- 1 PiV) 2 g- 1 iV o l] --{Q~ 1 No lj )VN+2N<j im oJ n = A^N, (34) 



3 



hij = -2Q- 1 iV ( j ij , (35) 

supplemented, of course, by Eq. (27). In other words, once Eq. (34) is 
solved for the trace- free part of the extrinsic-curvature tensor of S t , the 
induced metric on T, t is obtained from Eq. (35). This form of the coupled 
first-order ADM system of equations has possibly the merit of stressing the 
need to invert the operators L, Q and S to find a solution for given initial 
conditions. Approximate solutions, to the desired accuracy, will correspond 
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to the construction of approximate inverse operators L _1 , Q' 1 and S" 1 . For 
this purpose, it can be useful to re-express Eq. (34) in the form (I being 
the identity operator) 



I + -S~ 1 (L~ 1 VN) 2 Q~ 1 N 



{Q-^NoijyPN 



- S-\A ij N) = -2S- l (Na im ap). 



(36) 



The action of the inverse operator L 1 on VN is already given by Eq. (19), 

but the inverses Q~ l and 5 _1 make it necessary to develop an algorithm for 
their exact or approximate evaluation. More precisely, one has 



Q = L (i + ^L^NL^VNj , 
S = L(l - ^L^NL^TN^j , 



and hence 



Q- 



I + -L~ l NL~ lr PN 



I - -L^NL^TN 
3 



L~ 



U 



by virtue of the operator identity (AB) 1 = B 1 A 1 . 
view, two levels of approximation seem to emerge: 



(37) 
(38) 

(39) 
(40) 

From this point of 



(i) The degree of accuracy in the evaluation of Q 1 and S 1 . It is clear 
from Eqs. (34) and (35) that this reduces to finding (/ + pL~ 1 NL~ 1 VN) , 
where p = | or — |. A series expansion will be, in general, only of formal 
value. However, on using the norms defined in section 7.4 of Ref. [14], which 
relies on well known properties of Sobolev spaces, if the operator 



T p = pL^NL^VN 



(41) 



is a bounded operator on a Banach space X with norm ||T p || < 1, the 
inverse of I + Tp exists and the series J2ri=o(~ ^) n {T p ) n converges uniformly 
to (7 + Tp)^ 1 with respect to the norm on the set of bounded maps from X 
into X, and one can write 



(I + Tp)- l = J2(-l) n (T P ) n . 



(42) 



n=0 
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Equation (42) can be therefore used if L 1 and V are such that T p is a 
bounded operator with \\T p \\ < 1. 

(ii) The way in which the non-linear term OimO™ is dealt with. For example, 

one may regard the right-hand side of Eq. (36) as the "known term" , and 
hence consider the equation 



2 



lc-l 



Gadx 1 0dx J F S 

3 



{Q- l Na i:i )VN\ dx i <g> dx j 



- F~ 1 S- 1 {A ij N)dx i ® dx j = -2F- 1 S- 1 (Na im aJ n )dx i ® dx j , (43) 
where F _1 is the inverse of the operator 

F = I + -5 _1 (L _1 PAr) 2 Q _1 iV. (44) 

Equation (43) is an integral equation for (7jj, for which a perturbation ap- 
proach might be useful if one takes -2F- 1 5- 1 (iVa^o™) dx l <g> cix- 7 as the 

known term mentioned before. Note that we have resorted to the use of the 
tensor product dx l <g> dx^ because in the resulting equation (43) one has a 
well defined integration of a symmetric rank-two tensor field over the space- 
time manifold, here taken to be diffeomorphic to x R. Strictly, also Eq. 
(35) should be written as 

hijdx* <g> dx j = -2 (Q~ l Nvi^ dx i ® dx j , (45) 

and in all such equations only the symmetric part of the tensor product 
survives, because both hij and Gij are symmetric rank-two tensor fields. 

Note once more that it would be wrong to regard Eq. (27) as indepen- 
dent of the solution of Eqs. (28) and (29), because L^ 1 VN is only known 
at all times when the time evolution of hij and Kij has been determined for 
given initial conditions. Our equations (43) and (45), although not obviously 
more powerful than previous schemes, prepare the ground for an operator 
approach to the ADM equations, and hence might contribute to the un- 
derstanding of structural properties of general relativity. In particular, our 
way of writing the coupled system of non-linear evolution equations might 
lead to a better understanding of the interplay between elliptic operators on 
the spacelike hypersurfaces £t and hyperbolic equations on the space-time 
manifold. As far as we can see, this expectation is supported by the closed 
Friedmann-Lemaitre-Robertson-Walker model discussed before (where 
vanishes), and further examples of cosmological interest might be found, e.g. 
Bianchi IX models describing a closed but anisotropic universe. For open 
universes or asymptotically flat space-times, however, we are not aware of 
theorems that make it possible to expand the lapse as in Eq. (15). In such 
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cases, the counterpart of the elliptic theory on £ t advocated in our paper is 
therefore another open problem [15]. 
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